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Abstract 

We show, that the canonical invariant part of h corrections to the Gutzwiher 
trace formula and the Gutzwiller-Voros spectral determinant can be com- 
puted by the Bohr-Sommerfeld quantization rules, which usually apply for 
integrable systems. We argue that the information content of the classical 
action and stability can be used more effectively than in the usual treatment. 
We demonstrate the improvement of precision on the example of the three 
disk scattering system. 
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Gutzwiller trace formula for chaotic systems is often presented as the counterpart of 
the Bohr-Sommerfeld (BS) quantization of integrable systems Consequently, corrections 
of the trace formula proportional with powers of h are usually associated with quantum 
corrections [Q. In this letter we would like to show, that a part of h corrections is not 
connected to deep quantum effects and they can be calculated with some precision from 
purely semiclassical BS arguments. Then we propose a new trace formula and spectral 
determinant which is more precise than the usual trace formula although it uses only the 
linear stability and action as an input data, just like the original trace formula. 

First we would like to introduce periodic orbits from an unusual point of view. Chaotic 
and integrable systems on the level of periodic orbits are in fact not as different from each 
other as we might think. If we start orbits in the neighborhood of a periodic orbit and look 
at the picture on the Poincare section we can see a regular pattern. For stable periodic 
orbits the points form small ellipses around the center and for unstable orbits they form 
hyperbola. The motion close to a periodic orbits is regular in both cases. This is due to 
the fact, that we can linearize the Hamiltonian close to a periodic orbit, and linear systems 
are always integrable. Based on Poincare's idea, Arnold and coworkers have shown |Q, that 
the Hamiltonian close to a periodic orbit can be brought to a very practical form. One has 
to introduce new coordinates: one which is parallel with the orbit and others which 
are orthogonal. In the orthogonal directions we get linear equations. These equations with 
x\\ dependent rescaling can be transformed into normal coordinates so that we get tiny 
oscillators, or inverse oscillators, in the new coordinates with constant, frequencies. In the 
new coordinates, the Hamiltonian is 

1 "'"^ 1 

^ n=l ^ 

which is one possible normal form of the Hamiltonian in the neighborhood of a periodic 
orbit. The ± sign denotes, that for stable modes the oscillator potential is positive, while 
for an unstable mode it is negative. Since the eigenvalues of the monodromy or Jacobi 
stability matrix of a periodic orbit are invariant under the transformations we made sofar, 
the oscillator frequencies can be expressed for unstable modes with the Ljapunov exponent 
of the orbit 



cj„ = ln|Ap,„|/Tp, (2) 

where Ap „ is the expanding eigenvalue of the Jacobi matrix and Tp is the period of the orbit. 
Also, for stable directions the eigenvalues of the Jacobi matrix are connected with uj as 

Ap,n = e-'""^^ (3) 

The Hamiltonian (|l|) is integrable and can be semiclassically quantized by the BS rules. 
The result of the BS quantization for the oscillators gives the energy spectra 

En = (jn + fo^ Stable modes, (4) 
En = —ihun {in H — ^ ^^i Unstable modes. 
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where j„ = 0, 1, .... It is convenient to introduce the index s„ = 1 for stable and s„ = —i 
for unstable directions. The parallel mode can be quantized implicitly trough the classical 
action function : 

^ ^ p\ldx\\ = ^S\i{Em) = h (m + , (5) 

where fip is the Maslov index of the motion in the parallel direction. This latter condition 
can be rewritten in the equivalent form 

The eigenenergies of a semiclassically quantized periodic orbit are all the possible energies 

d-l 

E = Em + ^ En. (7) 

n=l 

This relation allows us to change in Em with the full energy minus the oscillator energies 
Em = E — J2n En- All the possible eigenenergies of the periodic orbit then are the zeroes of 
the expression 

A^^^{E)= Yl (1 - e*'^"^^"^"^'*"'^"^-'""^^/^^^/^"*'^^'"/^). (8) 
ji,---,jd-i 

Now, we show how one can derive the Gutzwiller trace formula from (|^). We have to 
expand the action around E to first order S\\{E + e) ~ Sp{E) + Tp{E)e, where Tp{E) and 
Sp{E) are the period and the action of the orbit, and we have to use the relations of a;-s and 
the eigenvalues of the Jacobi matrix, we get 

/ piSp{E)/h-iupn/2 \ 

^.W= n 1 - n n,v ' 

ji,...,jd-i V 

where i^p is the Maslov index of the orbit. Now, if we have many primitiv orbits and we 
would like to construct a function formally, which is zero whenever the energy coincides with 
the BS quantized energy of one of the periodic orbits, we have to take the product of these 
determinants: 

A{E) = llAp{E). (10) 

p 

This is exactly the Gutzwiller- Voros |1T0|J^ spectral determinant, which is the regularized 
semiclassical expression for the spectral determinant 

A{E) = det{E-H) 

of the Hamilton operator of the system. The logarithmic derivative of this quantity gives the 
trace of the Green's function and the oscillating part of the trace formula in semiclassical 
approximation 

J 1 pirSp{E)-iri/pT7/2 

TrGW, E) ^ -- ,og A(E) « - E T,iE)j^^^^—^, (11) 
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where the summation goes for the primitive periodic orbits and their repetitions, Jp is 
the monodromy or Jacobi matrix of the periodic orbit. The construction of the spectral 
determinant (|^) described above is not unique in the sense, that we can multiply each 
periodic orbit's spectral determinant with a smooth function which does not have zeroes 
and poles on the complex energy plane. However, the logarithmic derivative of such a 
smooth function will not contribute to the oscillating part of the trace formula and can be 
considered clS cl part of its smooth part. 

From this derivation we can see, that the Gutzwiller trace formula is recovered only if 
we linearize the action. This is a very bad approximation for low energies, where the ratio 
Tp{E) / Sp{E) ~ is large. In that regime we might expect that the spectral determinant 
(1) works better than the Gutzwiller- Voros formula @. 

That this is the case, we have checked on the example of the three disk scattering system 
13| , p^ at the standard parameters of disk separation 6 compared to the radius of the disks. 



This is a billiard system, where the parallel action is S{E) = kLp, where k = y/E is the 
wave number (with mass unit m = 1/2 ) and Lp is the geometrical length of the the orbit. 
There is only one unstable oscillator mode with oscillator energies 

^ ,, , ^ 2 A; log |Ar,L 

Epjik) = -ih + 1/2). 

Lp 

The new spectral determinant as a function of the wave number is 



A^^^{k) = — e*'^pV^^~^p^^/^"*''p'"/2). (12) 
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If we expand the exponent in h we can see that it gives corrections to the leading action and 
stability term. We can compare the first h correction with the exact h correction computed 
in Ref. |^ and Ref. The h corrected spectral determinant with the exact correction is 
defined as 



P j 



where C^^^^^^ the exact first h correction. The first correction coming from (p!2D is 



<r^'^ = (j + l/2)^(ln|Ap|)W 

On Fig. 1 we can see, that this accounts for about 80% of the exact correction. As we will 
see later, this is not the whole first h correction which comes form the BS quantization, it il- 
lustrates only that by using (0) we already take into account a whole series of h corrections. 
We also compared the exact quantum mechanical resonances ^ with those computed 
from the Gutzwiller- Voros spectral determinant and with the Gutzwiller- Voros spectral de- 
terminant with the first h correction |^ and plotted the results on the complex wave number 
plane. We can see on Fig. 2, that for large k the Gutzwiller- Voros spectral determinant, its 
h corrected version and the new BS type expressions approximate the resonances accurately, 
with a few percent error. However, the lowest resonances are approximated better by the 
new expression and the h corrected Gutzwiller- Voros determinant is even worse than the 
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uncorrected. This is because the Gutzwiller-Voros determinant and its corrected version are 
asymptotic series expanded in the powers of while the new formula approximates the 
eigenenergies of the individual periodic orbits for small values of k also correctly. We can 
conclude, that by using ([12D or in general we can considerably improve the Gutzwiller- 
Voros theory for low energies by using exactly the same input data (stability, action, Maslov 
index) in a more economic way. 

The semiclassical analysis outlined above can be done in a more general framework by 
the systematic application of the normal form theory of Birkhoff and Gustavson. This can 
be considered as a generalization of the normal form quantization of Swimm, Delos and 
Robnik [^] from equilibrium points to the periodic orbits. The full Hamiltonian expanded 
close to a periodic orbit can be written as the 'harmonic' plus an 'anharmonic' perturbation 

H{x\\,p\\,Xn,Pn) = Ho{x\l,p\\,Xn,Pn) + H a{x\\, Xn, Pn) , (14) 

where the anharmonic part can be written as a sum of homogeneous polynomials of x„ and 
Pn with x\\ dependent coefficients: 

oo 

if^(X||,X„,p„) = ^H^{x\\,Xn,Pn), (15) 
fc=3 

H\X\\,X^,P^)= hl^M)x^nPn"- (16) 

J2„ ln+mn=k 

We can carry out the canonical transformations introduced by Birkhoff and Gustavson 
pn],|8|JI|] and the Hamiltonian can be brought to normal form in the orthogonal directions. 
The X|| coordinate plays the role of a parameter. After the transformation up to order 
the Hamiltonian (|TBp is 

N 

^^(a;||,P||,Ti, ...rd_i) = Ho{xi\,pii,Ti,...,Td-i) + ^ f/^ (a;||, ri, rd_i) + h.o.t, (17) 

where is a j-th order homogeneous polynomial of r-s with x\\ dependent coefficients and 
Tn = ^{Pn ^ ^n^n) Hamiltoniau function of the original oscillator. This Hamiltonian 

truncated at order N is integrable, the nonintegrability is pushed to the higher order terms 
(h.o.t) . The orthogonal part can then be BS quantized by quantizing the individual oscil- 
lators, replacing r-s as we did in (|^). This leads to a one dimensional effective potential 
indexed by the quantum numbers ji, ...,jd-i 

H{x\\,p\\,ji, Jd_l) = -pI + U{x\\) + ^ hSnUJnUn + 1/2) + (18) 

^ n=l 

N 

+ J2 U^{x\\,hSiUJi{]i + 1/2), hS2UJ2{]2 + 1/2), hSd-lUJd-l{jd-l + 1/2)), 
k=2 

where jn can be any non-negative integer. The term with index k is proportional with 
due to the homogeneity of the polynomials. The parallel mode now can be BS quantized 
for any given set of j-s 
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Sp{E,ji, ...,jd-i) = j P\\dx\\ = (19) 

= j) dx\l^ 2 - ^ hSnUJnijn + 1/2) - t/e// (X|| , ji , jd_i)^ = 271% (m + ^ 



where Uejf contains all the x\\ dependent terms of the Hamiltonian. The spectral determinant 
becomes 



RSA 



(E) =Y\_ Yl (1 — e*'^*'*-^'-^^'"''-^''-^-'''^"*^'''^^^). (20) 
p iiv,id-i 



Here one can see, that the indices j, which were just auxiliary indices in in the Gutzwiller- 
Voros approach, now can be interpreted as "orthogonal quantum numbers". Formally, the 
integrability of the Hamiltonian is maintained in each level of the approximation and the 
Ruelle |jl2[ type zeta functions 



p 

corresponding to a given j configuration, can be interpreted as the factorization of the total 
spectral determinant according to subspaces with fixed quantum numbers. If we expand Sp 
in the exponent in the powers oi h Sp = J2k=o^''^k, we get corrections to the Gutzwiller- 
Voros spectral determinant in all powers of h. There is a very attracting feature of this 
semiclassical expansion, h in 5*^ shows up only in the combination hsnOJn{jn + 1/2). A 
term proportional with can only be a homogeneous expression of the oscillator energies 



Sn^n{jn + 1/2). We meutiou here Ref. [|l5l as a good example, where the superiority of the 
method outlined here can be demonstrated above the pure periodic orbit theory and the 
comparison of h expansion and BS quantization can be clearly studied. 

The h corrections derived here are doubly semiclassical, since they give semiclassical 
corrections to the semiclassical approximation. What can quantum mechanics add to this 
? Since quantum mechanics is not invariant for canonical transformations, the derived h 
corrections give only the leading behaviour of corrections and the exact corrections can be 
computed by other methods. The Birkhoff-Gustavson transformations should be replaced 
by quantum perturbation theory and semiclassical quantum numbers should be replaced 
by exact quantum numbers. This has been done in Ref. which we are going to publish 
elsewhere. We don't think, that the semiclassical determination of h corrections in higher 
orders is a very practical way to compute them, but in predicting the general structure of 
the corrections it helps us to understand their general behaviour and later probably to sum 
them up in order to get meaningful formulas also for low energies. 

The author thanks the discussions with E. Bogomolny, O. Bohigas, P. Cvitanovic, P. 
E. Rosenqvist, N. Whelan and A. Wirzba. This project has been financed by the EHCM 
PECO, OTKA F17166 and T17493. 
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FIGURES 



FIG. 1. Comparison between the exact h corrections and the BS approximation. We have 
taken all the periodic orbits p of the three disk system up to 9 bounces. Wc computed the first 
h correction for j = (Cpio)- Then we computed the quantity (In |Ap|)^/4iyp, which correction 
comes from the lowest order harmonic approximation of the BS quantization of the periodic orbit 
(Cp Q^^^^). We plotted here this value versus the exact h correction for each periodic orbit p. We 
can see, that for the 22 shortest orbits they are almost linearly correlated and about 80% of the 
exact correction comes from this effect. The rest is coming from the anharmonicity and from deep 
quantum effects. 

FIG. 2. Complex resonances of the 3 disk scattering system. The quantum calculation, the 
Gutzwiller-Voros approximation, the Gutzwiller-Voros approximation with one h correction added 
and the BS spectral determinant in harmonic approximation (BSH) . We used all the periodic orbits 
up to 9 bounces. For small Rek resonances the h correction breaks down, while our approximation 
is deviating the right way from the Gutzwiller-Voros result. 
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